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We construct explicit examples of microscopic models that stabilize a variety of fractionalized phases 
of strongly correlated systems in spatial dimension bigger than one, and in zero external magnetic 
field. These include models of charge fractionalization in boson-only systems, and various kinds of 
spin-charge separation in electronic systems. We determine the excitation spectrum and show the 
consistency with that expected from field theoretic descriptions of fractionalization. Our results are 
further substantiated by direct numerical calculation of the phase diagram of one of the models. 



I. INTRODUCTION 

Considerable theoretical effort has gone into under- 
standing the possibility of obtaining fractional quantum 
numbers for the excitations of strongly correlated sys- 
tems in two or more spatial dimensions, and in weak or 
zero magnetic ficlda^ Though much of the original inter- 
est arose in theoriesETLl of the high temperature supercon- 
ductors, ideas based on fractionalization have since been 
proposed to account for the properties of a numbet, of 
other poorly understood strongly correlated systemsoB. 
Field theoretic methods have enabled enormous-|Pp3|gress 
in obtaining a description of fractionalizationt3"E3El. A 
number of exotic fractionalized phases have been ac- 
cessed. The structure of the distinct possible gxcitations 
and the effective theory of their interactionsQ has been 
elucidated in some detail. A crucial feature is the pres- 
ence of discrete gapped vortex-like excitations - dubbed 
visons - apart from the excitations-.with fractional quan- 
tum numbers. It has become cleaiEd'Ej that fractionalized 
phases may be given a precise theoretical characterization 
through the notion of "topological order" — a concept elu- 
cidated clearly by WenllJ in work on the quantum Hall 
effect. 

Scepticism has been voiced in some quarters over these 
developments due to the almost complete lack of micro- 
scopic models that can be shown to display the phenom- 
ena mentioned above. Specifically, consider a model of 
a many-particle system with short-ranged interactions, 
and no special symmetries other than the global charge 
and/or spin conservation. Can fractionalization be shown 
to obtain in such a model? Apart from it's concep- 
tual value, answering this question will also help clar- 
ify the nature of the microscopic conditions that make 
it favourable for fractionalization to occur in a strongly 
correlated system. 

There has been limited (though important) progress in 
a direct answer to this question. Numerical studiesEl of 
a particular triangular lattice quantum spin-1/2 Heisen- 
berg magnet with ring exchange interactions provide evi- 
dence for a state with a spin gap and a four-fold degener- 
ate ground state on a torus as expected in a topologically 
ordered "spin-liquid" with fractionalized "spinan." excita- 
tions. In the context of quantum dimer modelai^, Moess- 



ner and Sondhilla have argued for a stable topologically 
ordered "liquid" phase on a triangular lattice. The stan- 
dard interpretation of the dimer model views the dimers 
as caricatures of singlet bonds formed between underly- 
ing Heisenberg spins on the lattice. From this point of 
view, the work of Ref. ^ provides supporting evidence, 
though not definitive proof, that models of Heisenberg 
antiferromagnets on triangular lattices do support a fraa- 
tionalized spin liquid phase. However, as is well-knownllZI, 
the quantum dimer model is exactly equivalent to a gauge 
theory — thus one may worry that establishing a topologi- 
cally ordered phase in the dimer model still docs not con- 
vince a sceptic that such phases can result in microscopic 
models with no special sviBmetries or a gauge structure. 

Recently, Balents et al£3 argued that a particular easy 
axis quantum spin-1/2 model on a Kagome lattice with 
short ranged (albeit complicated) interactions has a topo- 
logically ordered ground state with fractionalized exci- 
tations. This was done by reinterpreting it as a soluble 
point of the quantum dimer model on a triangular lattice 
but with three rather than one dimer emerging from each 
site, and following the same arguments as in Ref. (See 
also Ref. |l^ for a somewhat similar perspective.) How- 
ever, some features of this model such as the presence 
of two distinct visons, appear to be non-generic (from 
the point of view of the effective field theory of fraction- 
alized phases). This model also has an infinite number 
of local symmetries, and hence violates the requirement 
that fractionalization be demonstrated in models with no 
special symmetries. However, Balents et aZ.E3 made the 
important observation that perturbing the model slightly 
to get rid of the local symmetries will preserve the frac- 
tionalized phase. 

In this paper, inspired by these prior developments, 
we explicitly construct microscopic models that stabilize 
a wide variety of fractionalized phases. These include 
models for charge fractionalization in boson-only systems 
and various kinds of spin-charge separation in electronic 
systems. Our models involve only short-ranged interac- 
tions, and do not have any special symmetries other than 
global charge and/or spin conservation. We determine 
the excitation spectrum in the fractionalized phases and 
explicitly show the consistency with that expected from 
the effective field theories. Our results complete the an- 
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swer to the question of principle posed above, and will 
hopefully guide efforts to find materials and other even 
simpler models that realize these phases. 

We begin by illustrating our construction with a sim- 
ple Bose-Hubbard type model of a system of bosons (with 
charge ^f,) with short ranged interactions. The model has 
a global U{1) symmetry reflecting the conserved total bo- 
son number. We explicitly demonstrate the presence of 
two distinct Mott insulating phases in this model. In one, 
the excitations carry charges that are integer multiples of 
the underlying boson charge g^. In the other, there are 
excitations with boson number <7&/2, i.e. the bosons have 
fractionalized. This phase also has discrete Z2 vortices, 
the visons, which are gapped. Upon tuning a parameter 
in the model, it is possible to drive transitions from either 
Mott insulating phase to a superfluid phase. We further 
substantiate our arguments by performing a quantitative 
numerical calculation of the phase diagram of this model. 
The presence of topological order in one of the Mott insu- 
lating phases is detected numerically by the flux-trapping 
experiment discussed in Ref. We explicitly derive the 
effective field theory of the fractionalized phase and show 
that it is a theory of bosonic charge qb/2 fields coupled 
to a Z2 gauge field in it's deconfined phase. 

Next, we consider models of electrons coupled to su- 
perconducting phase fluctuations. These may be thought 
of as models of charge e electrons interacting with spin- 
less charge 2e bosonic Cooper pairs. We show how the 
boson only models above may be extended to include cou- 
pling to electrons to provide a realization of various spin- 
charge separated phases. These models therefore pro- 
vide explicit realizations of the routes explored in Ref. [ 
p^,p| for spin-charge separation. These spin-charge sep- 
arated phases have spin-0 charge e bosonic excitations 
(chargons or holons), spin- 1/2 charge neutral fermionic 
spinons, and a gapped vison. 

Of saecial interest are models that stabilize the nodal 
liquidtHlB (alias (1^2 _y2 RVB) phase — this has gapless 
fermionic nodal spinons, and has played an important 
role in theoppg-of the cuprate materials. While recent 
experimentf£ilO are not very encouraging on the possi- 
bility of fractionalization in the cuprates, it still is of the- 
oretical interest to demonstrate models that realize the 
nodal liquid phase. Another theoretically controversial 
possibility is that of ordered magnetic phases that nev- 
ertheless raije-|Spin-charge separated. This was also first 
discussedEJilj in the context of cuprate physics but is 
possibly relevant to a variety of other systems. We show 
how a model that stabilizes such ordered magnetic frac- 
tionalized phases may readily be obtained. This settles 
any doubts that may have been harbored on the possi- 
bility of such coexistence between magnetism and frac- 
tionalization. We then conclude with a brief discussion. 



II. FRACTIONALIZATION IN BOSON ONLY 
MODELS 



A. Model and general arguments 

Consider a system of bosons on the "face-centered" 
square lattice in 2D shown in Fig. |l| modeled by the 
Hamiltonian 

H = Hyj + -ffbond + -ffring + , (1) 

H^ = -W ^ {hl^rr' + h.C.) , 

-ffbond = -Jbond [(*rr')^(^^') + h.C. , 

{rr') 

Hrins = -i^ringXl(*l2*23*34*41 + h.C.) , 

□ 

r [rr') r 

Here bl = e'^'' are bosons residing on the corner sites 
of the lattice, and 'i'l^, = e^'^^^' are bosons on the bond- 
centered sites, which we identify by the end-points of 
the corresponding bond; n^, nf^, are the corresponding 
boson numbers, [6r,n^i] = iSrr', and similarly for ^'^r' 
and nfj.,. For technical convenience, we have chosen a 
rotor representation of the bosons (though this is not 
essential). The operator Nr is defined through 

Nr = 2n'^ + J2nt, . (2) 

r' Er 

The total boson number of the system is given by 

Ntot = lj2^-- (3) 

r 

The w-term is a boson hopping between the corner 
and the bond-centered sites, and r' G r sums over all 
such bonds emanating from r. The term i^ring is a ring 
exchange among four bond-centered sites belonging to 
the same square plaquette □, while the term Jbond is a 
similar ring-exchange-like boson interaction but among 
three sites associated with a given bond (rr'). The im- 
portance of ring exchange terms for promoting fraction- 
alization is stronglji suggested by the various field the- 
oretic descripi|ijms,EI and by previous studies of micro- 
scopic modelsJjH The Ub and are the usual on-site 
Hubbard terms. We have also included the Hubbard-JJ 
term for the boson number Nr. 

Despite the possibly unfamiliar form of the terms in the 
Hamiltonian, the following features are apparent. The 
model clearly has a global U{1) charge conservation sym- 
metry associated with global phase rotation of all the 
bosons. Note that if the 6-bosons are assigned charge 
Qb, then the ^f-bosons also have charge qb. There are 
no other special symmetries for general values of the pa- 
rameters. In particular, there are no local symmetries. 
Furthermore, all the interactions are short-ranged. We 
argue below that this model has a stable fractionalized 
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FIG. 1. Face centered square lattice on which our model 
Eq. dih is defined. 



insulating phase with charge ^ excitations and charge 
visons above a ground state with no conventional broken 
symmetries. 

Some gross features of the model can be guessed easily. 
At large the boson kinetic energy dominates and the 
system will be a superfluid. As w is reduced, there will 
be a transition to an insulating phase. The nature of this 
insulating phase depends on the other parameters in the 
model. In particular, the insulator will be fractionalized 
for Jbond: ^ringj large compared to ut,,u^. In the op- 
posite limit, a conventional Mott insulator with charge 
quantized in units of qi, will obtain. 

To establish these results, it is useful to consider the 
special limit w = 0; in this case, [Nr,H] = for every 
site r, and we can fix the value of iV^ for every r. Thus, 
in this limit, the model does have an infinite number of 
local symmetries. Later, we will move away from this 
special limit thereby destroying these local symmetries. 
For large U at w = 0, the ground state has Nr = 
everywhere. The w — model in the sector A^^ = for 
every r is readily understood as it can be regarded as the 
well-studiedE^ 3D compact ^(1) gauge theory coupled to 
a charge 2 scalar field. Indeed, divide the underlying 
square lattice of Fig. |l| into A and B sublattices. Let 

Or ^ Or = SrOr wlth 



= +1 if r e A , 
= -l if reB. 



(4) 
(5) 



To preserve the commutation relations, define the corre- 
sponding conjugate variables 



erUr 



(6) 



Similarly, let a^r' — (f'rr' iir Cz A,r' E B and a^r' = —(prr' 
iir G B,r' G A. Consider a as a vector field a^a = Or.r+aj 
with a = X, y, and perform the corresponding transfor- 
Era to the vector field Era conjugate to 



mation n^^, 
Urn- We have 
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FIG. 2. Phase diagram of the (2 -I- 1)D compact QED cou- 
pled to a charge 2 scalar. 

Nr = (A • E + 2n^) . (7) 
The Hamiltonian then becomes 

H = — 2 Jbond ^ COs{AaOr + 2ara) — '2K,-ing ^ COs( A X a) 



(8) 



while the constraint A^ = is simply the 'Gauss law' 

A • E + 2nl = . (9) 

As promised, H{w = 0) is the same Hamiltonian as for 
the (2 -I- 1)D compact QED coupled to a charge 2 scalar. 
This permits us to take over the classic results of Fradkin 
and Shenker on this model which determined the phase 
diagram to be of the form shown in Fig. ^ In the "con- 
fined" phase, all excitations carrying "gauge charge" are 
confined. In the "deconfined Higgs" phase, static exter- 
nal objects with gauge charge 1 are not confined. Fur- 
thermore, there is a stable gapped Z2 vortex (which we 
may identify with the vison) . A number of different per- 
spectives are available on these results. A useful physical 
one is to regard the deconfined Higgs phase as a "con- 
densate" of the charge-2 scalar. Naively, such a con- 
densate will have gapped vortices quantized in units of 
TT. However, due to the compactness of the gauge field, 
space-time monopoles are allowed in the theory. These 
correspond to events where the vorticity changes by 27r - 
consequently the vortices acquire a Z2 character. 

It is also clear that the deconfined Higgs phase has 
a topological order: e.g., the ground state is four- fold 
degenerate on a torus. These are simply obtained by 
threading no or one vison through the two holes of the 
torus. 

Consider now the excited states of the original Hamil- 
tonian H{w — 0) for large but finite U. Consider states 
such that Nrg = 1 at some site rp and A^^ = every- 
where else. Such a state can be regarded as a static gauge 
charge -1-1 at rg (assuming ro G A). In the confined phase 
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this sector costs an infinite energy in an infinite system. 
However, in the deconfined Higgs phase it costs only a 
finite energy. Remarkably, in the original boson model, 
such a state has a true electric charge of ^ (recall that 
Qtot = ^ X^r ^r)- Thus, in the deconfined Higgs phase, 
excitations with fractional quantum numbers for the true 
electric charge are allowed. In contrast, in the confined 
phase, finite energy excitations have gauge charge — 
this requires that J2rGA^r = J2reB'^r- Consequently, 
the excitations carry true electric charge that are integer 
multiples of qb, and hence are not fractionalized. 

These results on the w = Hamiltonian thus follow as 
a straightforward application of the standard Fradkin- 
Shenker analysis of the phase diagram of gauge theo- 
ries. However, they acquire even further importance here 
when we consider the Hamiltonian away from the w = 
limit when H no longer has an infinite number of local 
symmetries. Consider small w. This introduces fluctua- 
tions which mix states with different values of iV^ at the 
same site. However, for small w, these will not be capable 
of closing the gap to excitations about the ground state. 
Consider in particular the deconfined Higgs phase in the 
presence of a small w. The fractionally charged excita- 
tions are now allowed to hop from site to site and will 
acquire a kinetic energy of order w. However they will 
survive as meaningful excitations. The other indepen- 
dent excitation, namely the Z2 vortex, will also survive 
introduction of a small w. Thus, the original model has 
for w ^ but small a genuine fractionalized phase. (We 
can also add other more general boson hopping terms; 
clearly, the fractionalized phase will survive as long as 
these terms are weak.) In the subsequent subsections, we 
provide several direct confirmations of these arguments. 
In particular, we provide an explicit derivation of the ef- 
fective field theory of the fractionalized phase and show 
that it is a theory of charge ^ chargons coupled to a Z2 
gauge field in it's deconfined phase. This will also serve 
to make obvious our assertions on the properties of the 
model. 

We emphasize that despite the ease with which this 
result has been obtained, it has enormous significance. 
The Hamiltonian for w ^ has no special symmetries 
other than global charge conservation, and has only short 
ranged interactions. Nevertheless, it possesses a fraction- 
alized phase with charge ^ excitations and a gapped vi- 
son consistent with that expected from earlier field the- 
oretic descriptions of fractionalization. 

B. Numerical calculation of phase diagram 

In this subsection, we substantiate our results by a di- 
rect numerical calculation of the phase diagram of the 
model. To that end, it is useful to first consider a path 
integral representation of the model. The Euclidean ac- 
tion may be written 



S — e [Hw + -ffbond + -faring) 
r 

— ^ ^ COs(0rr-|-l — (^rr + 2Arr) 

rr 

— K'^ ^ COs(0rr',r + l — </'rr',r + \t + K't) 

{rr')r 

— iW^COaiXrr) ■ (10) 

rr 

Here, A^r G [0, 27r) is a phase variable living on the tem- 
poral links. To arrive at this form of the action, we 
first decoupled the Hubbard-?/ term in the path inte- 
gral and replaced all the Villain forms by cosines. The 
lattice spacing in the time direction is e, and the vari- 
ous couphngs are J'^ — l/(2Mf,e), K'^ = l/(2u^e), and 
W = l/(4C/e). 

The action represents a classical three-dimensional XY 
model with a global U{1) symmetry. As all the Boltz- 
mann weights are positive, we may analyse the phase 
diagram of the model using direct Monte Carlo simula- 
tions. To avoid unimportant complications, we will con- 
sider a particular choice of coupling constants where J"^ = 
2eJbo„d = J, K'' = 2eK,ine = K, and W = 2ew = W. 
Our choices of couplings J, K , and W are such that the 
resulting classical statistical mechanical system is rela- 
tively isotropic in space-timeEJ 

When W — the model is easily seen to reduce to 
the classical 3D compact QED coupled to a charge 2 
scalar. This has two phases neither of which has XY 
order (which implies insulating behaviour for the origi- 
nal quantum model) but which are topologically distinct. 
Turning on a small non-zero W does not induce XY or- 
der, but preserves the topological distinction between the 
two phases. Upon increasing W , there is eventually a 
transition to an XY ordered phase. Thus, we expect 
that a cut through the phase diagram in the K-W plane 
for large but finite J will look as in Fig. ^. 

We verify this expectation by direct simulations of the 
classical model Eq. ( |l0|) on cubic lattices of sizes up to 
12^^ with periodic boundary conditions. We use heat- 
bath local updates and run over 5, 000 Monte Carlo it- 
erations per each degree of freedom. We measure the 
XY order parameter in the original physical angles (e.g., 
M — Y^- e*^J ) and the superfluid stiffness ps associated 
with the direct boson hopping w [pa is defined in a stan- 
dard way — see, e.g., Refs. ^). Both quantities can be 
used to identify transitions into the superfluid phase. We 
also measure the specific heat of the classical system; this 
serves as an unbiased indication of the thermodynamic 
phase transitions and their order. From these studies 
performed at fixed moderately large J = 2.0, we obtain 
the phase diagram shown in Fig. |^, where we find three 
phases: an XY ordered phase {SC) and two distinct dis- 
ordered phases (T and I*). We should point out one 
detail about our scans through the parameter space: To 
perform an accurate study of the SC to X* transition, 
we always start from a fully ordered state inside the SC 
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FIG. 3. Phase diagram of the classical model Eq. (|lo|) at 
fixed J — 2.0. We label the phases using the language of the 
original quantum model Eq. (^): Superconductor SC is an 
XY ordered phase, while insulators T and I* are two mag- 
netically disordered but topologically distinct phases of the 
classical problem. 



phase. We found that if we start from a completely dis- 
ordered state in the I* phase, the system often traps a 
vison and subsequently a vortex when going into the SC 
phase (see also our discussion below), which significantly 
affects the measurements in our systems. 

We analyze the transitions using finite-size scaling. 
The J to SC transition for small K < 0.6 and the I* 
to SC transition exhibit 3D XY critical behavior; these 
are shown with open symbols in Fig. |^. For example, we 
can use the finite-size scaling relation, psL = g(L^^'^t), 
to locate the transitions and determine the correlation 
length exponent ly. Using this standard procedure, we 
also observe an important distinction between the two 
disordered phases: The universal value (psL)crit at the 
T* to SC transition is found to be one-fourth that at the 
I to SC transition, consistent with the charge fraction- 
alization in the X* phase. 

The I to SC transition for larger values of K approach- 
ing the T* phase, 0.7 < K < 1.0, seems to be first order; 
this is indicated with filled symbols and a heavy line in 
the same figure. Our evidence for this is the observed 
strong sharpening of the specific heat peak for the larger 
systems, with the maximum value growing very strongly 
with the system size. 

The 2 to X* transition (marked by crosses in Fig. ||) 
is most easily identified by observing the specific heat. 
This is a true thermodynamic transition, but is not ac- 
companied by any conventional ordering. The finite-size 
scaling of the specific heat peak is consistent with the 
3D Ising universality class. 

To illuminate the topological order in the I* phase, 
we perform flux trapping "experiments" as described in 
Ref. 20 (see also Ref. ^). We summarize these experi- 
ments in Fig. ^. The system is prepared deep in the SC 
phase with one vortex inside the annulus encircled by 
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FIG. 4. Flux trapping experiment at J = 2.0, K = 1.6. 
The system (of size 8^ ) is prepared with a single vortex in the 
annulus encircled by the Lx, deep in the SC phase (W = 0.7, 
cf. Fig. The vortex is detected by measuring the circu- 
lation Ix of the superfluid current (solid line). The system 
is cycled between the SC and X* phases. The "saw-tooth" 
dotted line is the MC time variation of W drawn so that the 
critical IVcrit coincides with the zero of I^- For a trapped vor- 
tex, the magnitude of Ix is set by the superfluid stiffness and 
is expected to be ~ ps-2n. The latter is shown with a dashed 
line for the first two cycles; the fact that the two quantities 
coincide indicates that the vortex remains trapped across the 
transition. 



the periodic L^', the physical angles 6 and (j) accumulate 
phase 27r going around the L^, and there is a superfluid 
current with circulation Ix = § j ■ dl « ps -^tt in this 
direction. As we decrease W towards the 2* phase, the 
vortex remains trapped all the way to the transition, and 
the magnitude of the superfluid current is set by ps- In 
the X* phase, the superfluid current is, of course, zero, 
but when we cycle the system back into the SC phase, 
the superfluid current reappears with full initial strength 
but with a random sign. For comparison, if we create 
a double vortex and perform a similar SC-T* cycle (not 
shown) on similar time scales, the double vortex "tunnels 
out" before we reach the X* phase and never reappears 
again. Similarly, a single or a double vortex created in 
the SC phase both disappear when we approach the X 
phase and never reappear again upon subsequent SC-X 
cycling. 

In terms of the effective degrees of freedom of the X* 
phase, the physical vortex is formed by a tt vortex in the 
chargon field and a vison. Bringing the system into the 
X* phase, the vison remains gapped and is trapped in the 
annulus. Cycling the system back into the SC phase, the 
vison binds a tt vortex in the chargon field, thus creating 
a physical vortex in the annulus but with a random sign. 
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C. Effective field theory 

We now provide a mapping of the model Hamiltonian 
to a Z2 gauge theory that wiU make obvious the results 
mentioned before. In addition, this yields an explicit 
derivation of the effective field theory for the fractional- 
ized phase. 

Consider the Hamiltonian Eq. ([|). To bring out the 
possibility of a fractionalized phase, define the operators 



g^j^j ^ gj</,,,^, through 



(11) 



Here = ±1 so that Bcr G [0, 27r). The field 6cr may be 
thought of as the "square root" of the operator hr and 
carries charge ^ and may be interpreted as a chargon 
operator. The field ^rr' on the other hand is charge 
neutral. 

Clearly, the boson number iV^ is conjugate to dcr and 
commutes with ^rr', while n^^^, is conjugate to ^rr' and 
commutes with Qrr'- 



[0cr,nt,]=O 



(12) 



We can now write the Hamiltonian in terms of 
{Ocr, Nr, cf>rr' rather than the original variables. 
However, to recover the original physical Hilbert space, 
we need to impose the constraint 



(13) 



this ensures that n^. = {Nr — J2r'er''^rr')/'^ integer 
(original 6-boson) number operator. 

Making this (exact) change of variables, we obtain for 
the parts of the Hamiltonian Eq. ([T]) 



Hhond 



r,r'^r 

-Jbond [(*rr')^ + ^-' 

{rr') 



(14) 
(15) 
(16) 



Note that the Ifbond term acts as an Ising anisotropy 
on the ^rr' field. Considerable simplification is possible 
in the limit of large Jbond and small Ub to which we now 
specialize. The potential "seen" by the phase (/)„' has 
two deep equivalent minima (j)rr' = or tt, which we label 
by cr^r' — e"^'-'' = ±1. The kinetic term {nf^,}'^ causes 
tunneling between the two wells. At each link, there 
are two low-energy states separated from all other states 
by a gap, leading to an effective 2-state system. In the 

cr^r' basis, we identify e^'^"rr' = a^^, since this operator 
translates (()„' by tt. Also, the kinetic term u^{n^^,)'^ 



is replaced by an effective transverse field ha^^, . In this 
large Jbond limit, the effective Hamiltonian becomes 

{rr') r 
-2ifring^Cri2CT23Cr|4Cr|i " ' ^^'^^ 

while the constraint Eq. ( |l3| ) is written as 

(-1)^" n = 1 (18) 

r' ^r 

at each site r. 

This effective model is precisely the quantum problem 
of chargons coupled to a fluctuating Z2 gauge field in 
two dimensions introduced and analyzed in Ref. |^. This 
model is known to have a phase diagram of the kind 
shown in Fig. ||. In particular, there is an insulating 
fractionalized phase where the chargon fields are decon- 
fined and there is a gapped vison (which occurs for large 
-ftTiing, and small w). 

D. Generalization to arbitrary commensurate filling 

Our results are readily generalized to arbitrary com- 
mensurate values of the total average number of bosons 
per unit cell. Consider a modification of the Hamiltonian 
where the U term is replaced by 



U{Nr - Nof 



(19) 



with A^o a constant. At such commensurate density (ra- 
tional values of Nq), insulating phases of the bosons will 
be possible. Again, in the limit of large Jbond, ^ring, 
and small this insulator will be fractionalized. 

Conventional (i.e. non-fractionalized) insulating states 
are of course possible in other limits. All of these cases 
are readily studied using the methods of the previous sub- 
section. Indeed, a non-zero A'o is trivially incorporated 
with no essential change leading to an effective Hamilto- 
nian Eq. (|l^) but with the modified U term above. As 
a special case of some interest, consider A'o pinLi The 
resulting model has previously been suggestedEIlEa as an 
effective model of frustrated easy-plane spin-1/2 quan- 
tum antiferromagnetsjp two dimensions. A recent study 
by Park and Sachdevn3 explicitly demonstrates the pres- 
ence of the expected two insulating phases: a bond den- 
sity wave crystal with confined excitations and a frac- 
tionalized phase. The fractionalized phase is more stable 
in this TVq = 1 case due to additional frustration coming 
from the Berry phase terms. 



III. MODELS FOR SPIN-CHARGE SEPARATION 

In this section, we generalize the models of Section || 
to construct models that display spin-charge separated 
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phases. We follow the route to spin-charge separation ex- 
plored in Ref. [ |l2|J| by considering models of electrons 
coupled to superconducting phase fluctuations. These 
may be thought of as models of spin- 1/2 charge e elec- 
trons interacting with spin-0 charge 2e Cooper pairs. As 
shown below, the independent excitations of the spin- 
charge separated phase are (i) a spin-0 charge e chargon, 
(ii) a spin-1/2 charge spinon, and (iii) a spinless charge 
neutral vortex — the vison. When either a chargon or 
spinon is taken all the way around a vison, the system 
acquires a phase of tt. This structure is exactly what 
is expected on the basis of the effective field theories of 
stable spin-charge separated phases. Indeed, as shown 
below, it is possible to provide an explicit derivation of 
the effective field theory as the correct description of our 
models in appropriate limits. 

In the models presented below, the spinous are 
fermions while the chargons are bosons. An important 
property of the spinous is that their number is not con- 
served. There are "pairing" terms in the Hamiltonian 
describing the spinon dynamics. Different spin-charge 
separated phases obtain based on the pairing symmetry 
of the spinous. Below we will discuss two different pairing 
symmetries as illustrative examples. 



A. Model for d-wave paired spinons 



Consider the following model: 



dwavc 



a,, 



{rr') 
(rr') 

(rr') 



h.c. 



(20) 
(21) 



K'( 



Cr-\Cr'l — CrlCr'\) + h.C. , (22) 



;7^(A^,- A^o) 



(23) 



Here, Cra represents the destruction operator for an elec- 
tron at site r and spin a. The electron is taken to have 
charge e. The operator '^rr' may, in this model, be con- 
sidered a Cooper pair living on the bonds of the lattice. 
In addition, these Cooper pairs on the bonds are coupled 
to other Cooper pair degrees of freedom br residing on 
the sites of the lattice. The corresponding boson-only 
terms H^, H^ond, and i?ring are the same as before. The 
operator A^^ is defined through 



Nr = 2nf, + 



(24) 



Clearly, the total charge Qtot = eJ2r'^r- The number 
A^o is a constant that sets the average charge per site. 
We take the "pairing amplitude" A^r' to have d^^-y^ 
symmetry. 

If U is large, the system will be in an insulating phase 
(for commensurate density) . The properties of this insu- 
lator depend on the values of the other parameters in the 



Hamiltonian. In particular, for large Jbond, ^ring, we ar- 
gue that the insulator will be spin-charge separated. The 
spinons are fermionic and have dx2-y2 pairing symmetry. 

We proceed as before and define the chargon field bcr 
and the neutral field ^rr' through Eqs. {Kw- It will 
also be extremely convenient to define a spinon field fra 
through 



Cra — ^crfra ■ 



(25) 



As before, the total charge associated with each site Nr 
is conjugate to the chargon phase 9cr and commutes with 

both ^!rr' and fra'- 



[Ocr,Nr] = I 



[*rr',A^r] = [/ra,Ar,] =0 . (26) 



As expected, the jra fields are formally charge neutral. 
Equations (O) also continue to hold. We further have 



[/ra,".^-'] = 



(27) 



and the equality c\.cfira — fra fra We may work with 
the set of variables {bcr, Nr, ^rr'jn^r'^ fra) instead of the 
original set {br,nr,'^rr','ntr'''^ra)- As with the boson- 
only models, this requires imposing a constraint on the 
Hilbert space which now takes the form 



(28) 



Continuing with the same steps as in Section ^ we find 
that in the large Jbond limit, the Hamiltonian reduces to 
the following: 



H = Ht+HA + H,i,[bc,a] , 
Ht^-tY, [flfr'blbcr' + h.c: 



(29) 
(30) 



{rr') 



Ha=J2 Kr'ifr^fr'l - /ri/r't) + h.c] . (31) 

(rr') 

Hchibc, 0"] is the same as before, Eq. (p7|). The constraint 
reduces to 



(32) 



We now argue that for Nq an integer, there is a sta- 
ble spin-charge separated phase. We first note that 
the Hamiltonian above describes a Z2 gauge theory of 
spinons and chargons coupled to the Z2 gauge field. As 
such, for large -faring 3> h, it's structure is almost iden- 
tical to the effective theory of a spin-charge separated 
phase of Ref. ^. The main difference is in the nature of 
the spinon hopping term (the term Ht) which seems to 
couple together the spinons and the chargons. However, 
this is readily seen to be an unimportant difference. 

First, consider the limit of small t,w (but t <^ w) at 
large repulsion U. In this limit, the chargons will lock 
into a Mott insulating phase (at integer A'o)- At t 
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1/7 = 0, the chargon number will be fixed at Nq per site. 
Going slightly away from this limit, we may treat both Ht 
and the chargon hopping term in perturbation theory to 
eliminate virtual charge fluctuations. The result will be 
an effective Hamiltonian describing the spinon and gauge 
degrees of freedom. To second order, the generated terms 
take the form 



E 

(rr') 



2U 



(33) 



with V — tflfr' + w<j^,^,. Expanding, we get two non- 
trivial terms: the first is simply spinon hopping coupled 
to the Z2 gauge field, while the second is a spinon four 
fermion interaction. The effective Hamiltonian then be- 
comes 

H = H,p^t+ ^^sp.int + fl"A + ^^IGT W] , (34) 

Hsp,t^-tspJ2<r'{fJjr'+h.C.) , (35) 

{rr') 

HspMt = -\J2[iPrfr'){fljr) + ifl,fr)if^fr')] ■ (36) 
(rr') 

Here the spinon hopping tgp — tw/U and the spinon in- 
teraction strength A = t^/{2U). Furthermore, in this 
large U limit, the constraint simply reduces to 



(-l)/.^/V+A'o ^ n <r- 



(37) 



r' £r 



As a function of -faring, this Hamiltonian undergoes a de- 
confinement transition. In particular, for large -faring, 
the fiuctuations of the gauge field may be ignored and 
the spinous are free to propagate. The nature of the 
spinon dispersion is easily found by considering the limit 
Kj-ing = 00. In this limit, we may set a^^, = -1-1 on every 
bond. The quadratic part of the spinon Hamiltonian is 
then formally the same as that describing non-interacting 
quasiparticles in a (1^2 _y2 superconductor, and therefore 
describes gapless nodal spinous. The spinon interaction 
is a formally irrelevant perturbation at this free spinon 
theory. As we are specifically in the limit that t ^ w, we 
have A ^ tgp — thus the interaction term may be safely 
ignored. Making Kj-i^g finite also only leads to irrelevant 
perturbations to the free spinon theory so that the long 
distance spin physics of the spin-charge separated phase 
is described by nodal fermionic spinous. 

The argument above considered the limit of large 
U, K^ing but small t, w. It is also instructive to consider 
the limit t^w ^ U . In this limit, the chargons are ex- 
pected to Bose condense leading to an ordinary €1^2 _y2 
superconductor. The long distance physics of this su- 
perconducting phase is readily captured by a continuum 
theory which keeps a continuum chargon phase field and 
the nodal spinous. The electron kinetic energy term Ht is 
then readily written as a spinon kinetic energy modified 
by the usual "Doppler shift" term coupling the gradient 



of the phase (the superflow) to a bilinear in the spinous. 
Vortices are permitted in this phase and have flux quan- 
tized in multiples of hc/2e. In the large i^ring limit, it 
is easy to see that the core energy of an hc/2e vortex 
will include a contribution proportional to ii'ring- On the 
other hand, the core energy of hc/e vortices does not di- 
verge as i^ring gocs to infinity. Now consider decreasing 
w to induce a transition to the insulator. At large idling, 
it is clear that this will occur due to proliferation of hc/e 
vortices rather than due to hc/2e vortices. Following the 
general arguments in Ref. , we will get a spin-charge 

separated phase. Note that, as argued in Ref. the 
Doppler shift term coupling the chargons and spinous is 
formally irrelevant, and one obtains a nodal liquid phase. 



B. Model for s-wave paired spinons 

It is straightforward to modify the model above to ob- 
tain one that stabilizes a spin-charge separated phase 
with s-wave paired fermionic spinons with a spin gap. 
We merely modify the pairing term above to 



Ha = ^E ^r^r^Crl + h.C. 



(38) 



Proceeding exactly as above, it is easily established that 
for large Jbond, C^, ^ring, such a spin-charge separated 
phase is indeed realized. 



C. Model for spin-charge separated magnetically 
ordered phases 

The effective field theories for spin-charge separated 
phases strongly suggest the theoretical possibility of spin- 
charge separation coexisting with magnetic long range 
order in a quantum phase. In this subsection, we show 
how the models above may be readily generalized to sta- 
bilize such phases. Consider a system consisting of two 
layers and a Hamiltonian of the form 



(39) 



Here H'-^\H'^^'> refer to parts of the Hamiltonian that 
depend only on the degrees of freedom residing in layers 
1 and 2 respectively. The interactions between the two 
layers are contained in the term H^^'^\ We assume that 
layer 1 consists of a square lattice of Heisenberg spins 
with 5 = 1/2 described by the Heisenberg antiferromag- 
netic model 



= J, 5: 

{rr') 



Sir ' ^Ir 



(40) 



We assume that layer 2 is described by the Hamiltonian 
Hdwavc in Eq. (|2^) above, and that the interaction be- 
tween the two layers is given by 
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J±_ } Si,. • (4^(TC2,.) 



(41) 



P. A. Fisher, Phys. Rev. B 62, 7850 



with ^ Ji. We assume that at Jj_ = 0, the layer 2 
is in it's spin-charge separated (and hence topologically 
ordered) phase. In this hmit, layer 1 will order antiferro- 
magnetically. Turning on a weak coupling J± will induce 
antiferromagnetic ordering in layer 2, but cannot destroy 
the vison gap. Consequently, the full model Hamiltonian 
will be in a phase that has magnetic long range order but 
nevertheless is spin-charge separated. 



IV. SUMMARY 

In this paper, we have discussed several concrete ex- 
amples of microscopic models in two spatial dimensions 
that display quantum phases with fractionalized excita- 
tions. These models possess no special symmetries other 
than those associated with global charge or spin conser- 
vation and also have only short ranged interactions, and 
thus confirm that fractionalization is a theoretically ac- 
ceptable possibility for strongly interacting many parti- 
cle systems in spatial dimension bigger than one. These 
models explicitly realize earlier field theoretic descrip- 
tions of fractionalization phenomena. 

A number of generalizations of our results are possible. 
Our models are easily generalized to arbitrary spatial di- 
mension, and provide concrete examples of fractionalized 
phases in any spatial dimension d > 1. For spin-charge 
separated phases of electronic systems, we have chosen 
to describe models with fermionic spinous and bosonic 
chargons. Following the ideas in Ref. ^ these are read- 
ily modified to construct spin-charge separated phases 
with fermionic chargons and bosonic spinous (at least 
with easy plane spin anisotropy). An additional upshot 
of our results is the couRttuction of topologically ordered 
classical 3D XY modelscJ. 

Finally, we mention that quantum nliases with topo- 
logical order have also been suggestedtHl to be suitable 
states of interest to quantum computation. The topolog- 
ical structure naturally protects the systeimj^om deco- 
herence. This very preliminary applicatiorJlS'LlI may also 
benefit from the results in this paper. 

We thank Leon Balents, Matthew Fisher, and Subir 
Sachdev for useful conversations. This work was sup- 
ported by the MRSEC program of the National Science 
Foundation under grant DMR-9808941. 
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